Mixing of magnetic and phononic excitations in incommensurate Spin-Peierls systems 
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We analyze the excitation spectra of a spin-phonon cou- 
pled chain in the presence of a soliton. This is taken as a 
microscopic model of a Spin-Peierls material placed in a high 
magnetic field. We show, by using a semiclassical approx- 
imation in the bosonized representation of the spins that a 
trapped magnetic state obtained in the adiabatic approxima- 
tion is destroyed by dynamical phonons. Low energy states 
are phonons trapped by the soliton. When the magnetic gap 
is smaller than the phonon frequencies the only low energy 
state is a mixed magneto-phonon state with the energy of the 
gap. We emphasize that our results are relevant for the Ra- 
man spectra of the inorganic Spin-Peierls material CuGeOa . 



The discovery in 1993 by Hase et al H of the first in- 
organic Spin-Peierls compound CuGeOs has opened the 
possibility of study the physics of this collective phenom- 
ena in a deep way. Several experimental proves have 
given an exhaustive information about the excitation 
spectra of this system and its evolution with an applied 
magnetic field. The effect of non-magnetic impurities has 
been investigated also. 

Theoretical studies have focused on simplified mag- 
netic model. The excitation spectra in the low tem- 
perature phase have been analyzed using a dimerized 
and frustrated Heisenberg chain as a minimal model for 
this material [||-|[. The logic underlying these studies 
are: the competition between magnetic and elastic ener- 
gies resolves in the low temperature phase in the dimer- 
ization of the lattice. Once this process takes place, 
phononic and magnetic excitation completely decouple 
and the magnetic excitations are the same as the chem- 
ical dimerized system. This point of view is based on 
an adiabatic approximation supposing that the energy 
scale of the magnetic process are high enough respect to 
the phononic ones. As it has been recently emphasized 
[pi, this relation is not fulfilled for CuGeOa where the 
phonons relevant for the dimerization process are about 
one order of magnitude more energetic than the mag- 
netic gap. The adiabatic approximation is questionable 
for this system. An antiadiabatic approach has been de- 
veloped. The frustrated interaction arises, in this con- 



text, from the integration of the in-chain phonons and 
the explicit dimerization from the interchain interaction 
treated in a mean field approximation |7|. The same 
frustrated-dimerized Hamiltonian is therefore obtained 
but with a reinterpretation of the parameters. What is 
clearly missed in these studies is a general understanding 
on how spin and phonons mix as elementary excitation 
and how the spectra of Spin-Peierls systems is built as a 
result of this mixing. Some recent numerical results have 
partially addressed this question [|j. 

In this paper, we analyze the excitation spectra of one- 
dimensional spin-phonon system by semiclassical tech- 
niques on the bosonized representation of the spins sub- 
system. We focus on the properties of this system in a 
high magnetic field. In the dimerized phase the system 
is in a singlet ground state. Coupling with a magnetic 
field is not effective up to a critical field value where 
some dimer break and their spins are freely from their 
singlet. The lattice also relaxes in this process forming 
the so called soliton lattice. This relaxation is a first 
indication that spin and phonons do not act as indepen- 
dent excitation. It has been experimentally shown by RX 
measurements that the lattice becomes incommensurate 
following a soliton pattern P]. The magnetic profile of 
the soliton has been analyzed by NMR measurement [Q . 

More recently, optical proves have shown the spec- 
tral signature of the incommensurate phase. In the uni- 
form dimerized phase a low energy resonance appears at 
30 cm~^ which is a smaller energy than two times the 
magnon gap [nl|. This peak has been adjudicated to a 
bound state of two triplet M. When the incommensu- 
rate phase appears the spectral weight is transfcrcd from 
this state to a lower energy peak at the position of the 
magnon gap O] . Direct magnon process could not be seen 
in the optical response due to spin conservation. Soliton 
assisted one magnon excitation was recently been pro- 
posed as the origin of this peak ^ in similarity with the 
situation in the presence of non-magnetic impurities [[l2| . 

This approach is based on the image of the soliton as an 
isolated spin in the externally dimerized chain so that an 
adiabatic approximation is supposed. We will show that 
this state in fact disappears when the dynamics of the 
phonons is included . We will also show that a trapped 



magneto-phonon state with the energy of the gap appears 
in the antiadabatic regime so explaining the optical data. 
Let us proceed more formally. Our starting Hamilto- 
nian for an Spin-Peierls compound is: 
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Si are spin-1/2 operators of the i ion and a the magneto- 
elastic coupling . Hph represents our simplified model for 
the phonons. It contains a scalar coordinate Ui and its 
conjugate momentum Pi which are supposed to be the 
relevant ion coordinates for the dimerization process. As 
we will only retain the phonons relevant for SP transi- 
tion, the specific dispersion of the phonons is not im- 
portant in our approach. The tridimensional character 
of the phonons field is essential to account for the finite 
Spin-Peierls temperature and the excitations in the low 
temperature dimerized phase |l3[. We will discuss later 
its effect on the incommensurate phase. 

The low energy spectrum could be studied by bosoniza- 
tion. The spin variables are approximately represented 
by the bosonic field (f>{x) and its conjugated momentum 
n(a;) |QJ|. In addition, we retain only the phonon modes 
producing a smooth deviation of the dimerized pattern. 
So, we make the replacement (— l)*Wi — > u{x). The low 
energy Hamiltonian becomes: 
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where a is the lattice constant of the original chain and 
Uq a short range cutoff introduced in the bosonization 
procedure. Vg is the spin wave velocity and 77 is related 
with the exponent of the correlation functions. For the 
isotropic Heisenberg model with nn interaction we have 
^, = ^ and ?7 = 2. 

We use a semiclassical approach in the form of a self 
consistent harmonic approximation(SIIA) as it has been 
originally proposed by Nakano and Fukuyama (NF) [g5| . 
In this approach the boson field is split in a classi- 
cal component (j)ci and the quantum fluctuation (f). The 
last term of (Q) is developed up to second order around 
(j)ci and then treated self-consistently by the following re- 
placement: 
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(02) is the ground state expectation value. 

In their original work NF have fixed the displacement 
field u to its equilibrium classical value {ud), so that an 
adiabatic approximation was assumed. Let us summarize 
the main results of this study and its consequence for the 
spectra: 

• The classical equations for and u{x) have a soli- 
ton solution of the form: sin(0c;) = tanh(a;/^) 
Uci{x) = uotanh(x/^) . ^, the soliton width, 
is Ws/A and A is the gap over the homogeneous 
dimerized state (the 'magnon' gap) and uq are the 
displacements in this homogeneous dimerized state. 
The soliton carries Sz = f /2 spin. In the presence 
of an external magnetic field these solitons will con- 
dense in the ground state and accommodate in a 
soliton lattice structure. 

• The creation energy of the soliton as well as the ex- 
citations over this solitonic vacuum could be eval- 
uated once the eigenvalues of the fluctuation op- 
erator are known. This eigenvalue problem cor- 
responds to a one of a Schroedinger-like equation 
which reads as: 



{-V? +tanh^(2:)}V'A = 



n{\f 



V'A(i) 
x 



(6) 



f2(A) are the frequency oscillations in presence of 
the soliton. Eq. 



has one bound state at ^ 



75-1 



0.786 and a continuum started at Q. 



. The excitation spectra of the theory is spanned 
by the following state: 

• The ground state of the quantum soliton builded 
around 0c/. Their creation energy {Es) is given 
by its classical energy plus the difference between 
the sum of the zero point energies of the oscillators 
(rj(A)) and the ones in absence of the soliton. This 
creation energy measures the critical field for the 
commensurate- incommensurate transition [n5| . 

• The excited state of the soliton with energy E^ = 
Es + ^b- It is related with the bound state of Eq. 

• Labeling by q the continuum of level of Eq. ||, 
they are ^{q) — Vs\/q'^ + I/^^ This is just the ki- 
netic energy of a particle in the soliton free sector. 
In terms of the original spin chain this is the low 
energy dispersion of the Sz = component of a 



magnon. Therefore this state corresponds to the 
scattering of a magnon in presence of the soliton. 
Moreover, when one of the continuum modes is ex- 
cited once, we get a two particle magnon-sohton 
state. 

The two last adjudication of states have not been done 
in the original analysis of the NF work. As our formalism 
breaks SU(2) symmetry at intermediate stage, the total 
spin of each state is not directly accessible and should be 
carefully reanalyzed. We have previously stated that the 
states below the continuum correspond to the emission of 
a magnon in the presence of the soliton. Therefore they 
have total spin-3/2. The bound state has the same total 
spin as the soliton. Then, it is an optical active mode 
and it is the analogue in our formalism of the so called 
soliton assisted magnon process founded in the strong 
dimerization limit |12| . 

Moreover, the previous predicted spectra could be 
tested by numerical exact diagonalization on finite chain. 
To this end, we solve iteratively the adiabatic equations 
for Ui arising from Hamiltonian Eq. |l|. The ground state 
of the spin system was recalculated at each iteration step 
by Lanczos diagonalization. We considered chains of odd 
number of site up to L = 23 sites with periodic boundary 
conditions and look for the equilibrium positions in the 
subspace Sz — 1/2. The numerical details of the method 
was given in [ pTlJl^ . 

Once the ionic coordinate converge we diagonalize the 
spin Hamiltonian, with fixed Ui, to obtain the low lying 
excitation in the Sz = 1/2 and Sz = 3/2 subspaces. The 
results are shown in Fig.H as a function of the 1/L^. Pa- 
rameter (^M = \ has been chosen in order to have a 
thin soliton and to reduce finite size effect. The bound 
state predicted by the bosonized theory is clearly seen. 
The higher energy states collapse in a continuum of total 
spin-3/2 in the thermodynamical limit. By linear ex- 
trapolation, we found for L —>■ oo the ratio 0.84 between 
the bound state and the border of the continuum. This 
value compares well with our previous prediction 0.786. 
We conclude that the method gives at least a qualitative 
features of the low energy spectra. The states appearing 
in the border of the continuum in Fig. (|lj) are an artifact 
of the strongly thin soliton we are considering. We have 
check that for bigger values of -t^jy^ ~ j these states in 
fact disappear. 
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FIG. 1. The excitation spectrum calculated by Lanczos 
diagonalization in the presence of an adiabatic soliton. Circles 
are for states of 5*2 = | and crosses for Sz = ^- The zero of 
the energies is chosen at the soliton ground state. Energies 
are in unit of J. 

We discuss in the following two effects that could 
change our estimated value of the bound state from the 
SHA method: 

1. In going from |^ to Q, we have neglected as NF in 
their original work a term proportional to cos(20) 
in the bosonized theory. Even though this oper- 
ator plays an important role in the case of Ising 
anisotropy, it is a marginally irrelevant operator 
in the Heisenberg case and therefore it should 
not change the qualitative features previously dis- 
cussed. Moreover, NF have studied in a further 
contribution ||l^ the effect of this term. Unfortu- 
nately, the estimation of the bound state value is 
strongly dependent on the value of the cuttoff in 
the bosonization procedure (oq) in this situation. 
If we assume as NF that Trag ~ a, the bound state 

becomes VlhA = \ hQ- ^ 0.93. We can conclude as 



a general fact that the presence of this term tend 
to increase the bound state value. 



2. It has been recently noted [^J20| that a full self- 
consistent treatment (where < (fp' > m Eq. is not 
fixed to its value in the ground state) changes the 
form of the soliton. Moreover, it has been shown in 
Ref.( iQ) that this treatment produces a smoother 
soliton than the one considered here. Therefore, we 
conclude that this treatment will reduce the bound 
state value. 

In the following we analyze the evolution of the spectra 
when non-adiabatic effects are included. For simplicity 
we do not introduce the two additional effects previously 
discussed. These effects will not affect the prediction of 
the trapped magneto-phonon state we will find. 



Inclusion of non-adiabatic effects. To go beyond this 
static adiabatic approximation in this semi-classical cal- 
culation it is necessary to include fluctuations in the dis- 
placement field u{x). Therefore, we split u{x) as Ud +u 
and replace it in Eq.(Q). The classical equations are the 
same as before. The fluctuation operator is now a 2 x 2 
differential operator with component over u and (j). The 
eigenvalue problem is now given by: 
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■01 and "02 are the component of the fluctuation eigenvec- 
tor over (j) and u respectively, w = y -^ is the phonon 
frequency at q = tt The eigenfrequency are as previously 
given by fl. We obtain from the second equation: 



■02 



uoa;^sech(x) 



01 



(8) 



Replacing in the flrst equation we have a kind of 
Schroedinger equation where the potential depends on 
the energy. This equation reads: 
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It is worthy to note that, the only relevant parameter is 
■g-.Then, we can follow the evolution of the spectra from 
the adiabatic to the antiadiabatic limit by moving this 
parameter. 

We have numerically solve Eq. (g). In Fig. ||, we show 
the evolution of the excitation spectra as a function of -^ . 
The soliton creation energy has been taken as the zero of 
the energy. 




FIG. 2. The excitation spectrum as a function of the 
phonon frequency. Both energy scales are given in units of 
the magnon gap (A) 



The main results are: 

• There is a zero energy solution which was not 
present in the adiabatic calculation. This mode 
arises in the complete translation invariance of the 
theory. It is nothing but the translational mode 
of the soliton. In a realistic situation solitons are 
not isolated, they form a soliton lattice. This zero 
mode will acquire some dispersion giving rise to the 
so called phason modes of the soliton lattice. If the 
magnetic field is not much higher than the critical 
one the solitons will be rather separate and we can 
neglect their interference. 

In the direction perpendicular to the magnetic 
chains, solitons will accommodate in a domain wall, 
i.e. they form an array of parallel solitons. Their 
coherency is assured by the interchain elastic cou- 
pling. The zero mode as well as all the states of 
the spectra will have a transversal dispersion, they 
will depend on the momentum perpendicular to the 
chains. The spectrum shown in Fig. (0) gives the 
q=0 excitations. We will analyze their dispersive 
character in a forthcoming work. 

• The upper branch corresponds to the excited state 
of the soliton previously found in the adiabatic ap- 
proximation. Its energy increases with w up to a 

2 1 

critical value ^ = ^ where this state disapear. For 
this w and $1 = A, the energy dependent potential 
vanishes. The lattice soliton produces an harmonic 
potential over the magnetic soliton where it oscil- 
lates. In the adiabatic approximation the lattice 
soliton is frozen and the magnetic one oscillates. 
For finite u, the lattice soliton can also move and 
vibrate with respect to the magnetic deformation. 
Finally, at the critical value of w there is no more 
a localized vibration and the soliton becomes an 
unique entity both for the lattice and for the spins. 
This explains the disparition of the trapped state 
at a critical phonon frequency. 

• The straight line of slope one in Fig.y corresponds 
to fi = u. The eigenstate with this energy corre- 
sponds to the excitation of a phonon in the pres- 
ence of the soliton. The eigenvalues near above 
are satellite phonons corresponding to phonons 
trapped by the soliton (i.e. particle vibrations near 
the soliton). They have only spectral weight on 
the phonons as it can be seen putting in Eq.(g) 
n ^ to. For increasing phonon frequency, most of 
these eigenvalues lose but the lowest survives ac- 
quiring the energy of the magnetic gap. This state 
is a mixing of a magnetic and a phononic excitation 
as we show in the following. 

• Four Peierls active phonons have been identified for 
CuGeOa pi]]. They have frequencies of 3.12, 6.53, 



11.1 and 24.6 THz i.e. the smaller one is about 
150K. The spin gap is rather small for this system 
A — 24:K. As we previously stated the real pa- 
rameter regime where this material lives is w > A 
and the only survival bound state has the energy of 
the gap. This is a singlet state which we associate 
with the low energy peak seen in Raman scatter- 
ing in the incommensurate phase 0. The ipi and 
rp2 components of the corresponding eingenfunction 
are different from zero for this state. 

Their mixed character could be advocated by analyzing 
separately the magnetic and phononic spectral response. 
The magnetic Raman operator of a dimerized chain is 
given by 7^™g = J2ii^ + 7mg(-l)'Si • Si+i) ^ j„ig is a 
microscopic parameter. In our bosonized formalism the 
most relevant contribution comes from the staggered part 
and it is given by ^^^^^^sm(0). By using the SHA given 

by Eq. (ra) , retaining the term linear in (j> and developing 
in the basis of the eigenvalues of the fluctuation operator 
the magnetic spectral weight (M„) of an state of energy 
fin is given by: 



71 y- _ ''»g 



AK Qr 



Jaa 



ip2dx 



(10) 



where we have used Eq. (H) to write ■01 as a function of 
^p2■ Phonons contribution to Raman scattering is propor- 
tional to the square of the transition elements of the nor- 
mal coordinates TZph = 7ph J ~^u{x), the phonons spec- 
tral weight of the n-state is given by: 
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The relation between the magnetic and phononic con- 
tribution to the n-th peak of the spectra is: 
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We take the first two factors as a measure of the spin- 
phonon coupling , we fix this parameters when moving w. 
Eq. (|2|) shows that the magnetic component of a given 
Raman peak increases when its position shifts from a 
phonon frequency. The lower energy eigenvalue of Fig(0) 
starts being a trapped phonons for -^ < 1 and transmutes 
in a magneto-phonon in the opposite limit where the real 
material lives. Its existence is a consequence of the quasi- 
onedimensional character of our system because a one- 
dimensional well has always one bound state. 

In summary, we have shown that nonadiabatic effects 
are relevant for incommensurate Spin-Peierls system. In 
the antiadiadatic limit applicable to CuGeOs we show 
the apparition of a trapped mixed state with the energy 
of the magnon gap. It gives an alternative explanation of 
the Raman induced peak found in the incommensurate 
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